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The  Statistical  Relations  Among  Q,  Kp,  and 
the  Global  Weather  Central  -K  Indices 


1.  INTRODUCTION 


In  the  geographical  region  covered  by  the  Over -the -Horizon  (OTH)  Experimen¬ 
tal  Radar  System,  the  auroral  ionosphere  is  one  of  the  most  important  environmen¬ 
tal  factors  affecting  the  performance  of  the  radar.  It  is  known  that  the  high  latitude 
ionosphere  and  the  auroral  oval  are  best  ordered  by  the  high  latitude  magnetic  activ 
ity  index  Q.  The  diameter  of  the  auroral  oval  and  its  location  are  best  expressed 
in  terms  of  the  Q  index.  On  the  other  hand,  high  latitude  ionospheric  phenomena 
have  been  expressed  as  a  function  of  global  geomagnetic  activity  using  the  Kp  index. 
Therefore,  to  provide  the  future  Operational  Radar  System  with  specifications  of 
the  auroral  and  sub-auroral  (F-trough)  ionosphere,  it  is  important  to  have  Q  and  Kp 
indices,  or  other  equivalent  parameters,  available  in  real  time. 

Ideally,  the  Q  index  is  determined  by  a  measurement  on  the  midnight  magnetic 
meridian  at  a  geomagnetic  latitude  A  >:  58°.  At  present,  the  Q  index  is  measured 
at  only  one  station,  Sodankyla,  Finland. 

The  Q  index  measured  at  a  single  station  shows  a  typical  diurnal  variation1  due 
to  the  movement  of  the  station  with  respect  to  the  auroral  oval  (explained  in  the 
analysis  section).  This  systematic  diurnal  variation  prevents  the  use  of  a 

(Received  for  publication  11  January  1982) 

1.  Dandekar,  B.S.  (1979)  Magnetic  Disturbance  Statistics  from  a  Single  Station  Q 


single -station  Q  index  to  specify  the  position  of  the  auroral  oval  at  all  Universal 

Times.  The  planetary  index  of  magnetic  activity  Kp  can  be  used  to  specify  the 

position  of  the  auroral  oval.  However,  a  delay  averaging  2  months  occurs  before 

the  Kp  index  becomes  available,  .which  prevents  the  use  of  the  Kp  index  for  real 

time  prediction  or  specification  of  the  auroral  oval. 

Because  no  index  of  global  magnetic  activity  is  available  on  a  real-time  basis, 

the  Air  Force  Global  Weather  Central  (AFGWC)  relies  on  a  determination  of  a 

3 -hourly  K  index  computed  from  the  magnetic  observations  from  five  observatories. 
2 

A  preliminary  analysis  based  on  a  small  sample  of  data  from  4  months  (June- 
October  1980),  showed  the  feasibility  of  using  the  AFGWC  K  index  in  deriving  Q 
(or  Kp)  required  for  the  specification  of  the  auroral  ionosphere.  This  report  pre¬ 
sents  the  results  from  a  much  wider  data  base  of  39  months.  For  convenience,  the 
AFGWC  K  index  will  be  referred  to  simply  as  the  K  index. 

The  study  showed  that  the  K  index  tracks  the  Kp  index  fairly  well.  However, 
the  K  index  has  some  shortcomings:  For  low  values  of  Kp  the  K  index  is  l/3  unit 
larger  than  the  Kp  index;  for  high  values  of  Kp,  the  K  index  is  l/3  unit  smaller  than 
the  Kp  index;  and  since  small  values  of  K  and  Kp  are  more  common,  the  K  index 
is  l/3  unit  larger  than  the  Kp  index  for  most  values  measured.  In  addition,  the 
K  index  exhibits  some  diurnal  variation.  Empirical,  time  dependent  factors  have 
been  determined  to  correct  K  so  that  it  more  closely  corresponds  to  Kp.  Also, 
empirical  relations  between  K,  Kp,  and  Q  have  been  determined  for  obtaining  the 
Q  index  from  K. 

2.  DATA  BASE 

The  AFGWC  K-index  data  for  March  1978-May  1981  were  used  in  this  study. 
The  Kp  index  and  the  Q  index  from  Sodankyla  were  also  used  in  the  analysis.  The 
latter  data  sets  are  from  the  World  Data  Center  A,  Boulder,  Colorado. 

Instead  of  the  conventional  quasilogarithmic  K  index,  AFGWC  provided  the 

3 

linear  ap  index.  These  data  were  converted  to  the  AFGWC  K  indices.  In  the  fol¬ 
lowing  analysis  it  is  assumed  that  all  these  data  have  been  properly  scaled  and  need 
no  correction  for  observational  procedures. 

During  the  period  of  39  months,  AFGWC  implemented  several  changes  in  the 
K-index  network  operation.  Therefore  the  data  were  initially  divided  into  four 
corresponding  separate  sets  for  analysis.  Since  no  significant  differences  were 

2.  Dandekar,  B.S.  (1979)  Study  of  the  Equatorward  Edge  of  the  Auroral  Oval  from 

Satellite  Observations,  AFGL-TR  -79-0010,  AD  A072  997. 

3.  Prochaska,  R.D. ,  Capt. ,  Geomagnetic  Index  Calculation  and  Use  at  AFGWC, 

AFGWC/Tn-80/002 . 
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found  among  the  separate  data  sets,  all  the  data  were  treated  as  a  single  continuous 
set  in  the  following  analysis. 


3.  ANALYSIS 

The  averaged  K  determined  for  all  data  taken  for  a  particular  Kp  level  should 
correspond  exactly  to  that  particular  K.  A  measure  of  the  quality  of  the  K  index  is 
how  well  K  matches  Kp  for  a  selected  Kp  level  and  the  standard  deviation  a  of  the 
respective  K  distribution.  A  large  a  indicates  the  data  base  may  be  too  small  for 
K  to  systematically  produce  a  reliable  Kp  estimation.  A  systematic  difference  be  ¬ 
tween  average  K  and  Kp  would  suggest  a  weighting  error  of  the  input  (station)  K 
values,  or  neglect  of  possible  diurnal  variations  emphasized  by  the  small  longitude 
sector  (see  Figure  2,  which  is  shown  later)  providing  inputs  to  K,  to  systematically 
produce  a  reliable  Kp  estimation. 

All  this  can  be  tested  by  the  mass  plot  of  the  data.  In  practice,  the  fixed  inter¬ 
vals  for  the  indices  combined  with  the  large  number  of  data  points  results  in  a 
thick  superposition  of  about  150  separate  points.  This  limitation  can  be  bypassed 
by  studying  the  observed  frequency  distributions  of  the  observations  in  the  corre¬ 
sponding  intervals  of  K  and  Kp.  This  distribution  is  presented  in  Table  1.  In 
Table  1  the  observed  frequencies  in  the  Kp  intervals  are  along  the  columns  and  the 
frequencies  in  the  K  intervals  are  along  the  rows.  The  sum  of  the  occurrence  fre¬ 
quencies  along  the  row  is  in  the  last  column,  and  the  sum  down  the  column  is  in 
the  bottom  row. 

In  Table  1,  the  lowest  value  K  =  0Q  (last  column,  row  1,  summation  frequency 
=  0)  is  absent.  The  corresponding  lowest  value  Kp  =  0Q  is  seen  139  times  (in  the 
first  column  and  last  row).  The  highest  value  Kp  =  8+  corresponds  to  the  highest 
value  K  =  9  .  The  median  range  (50  percent  of  the  observations)  of  Kp  =  2  cor¬ 
responds  to  K  =  2+.  The  highest  observed  frequency  of  Kp  is  between  2+  and  3  , 
whereas  the  corresponding  range  of  K  is  2_  to  2+.  The  results  show  that  K  is  a 
reasonably  good,  but  not  identical,  replacement  foi’  the  Kp  index. 

To  determine  the  correlation  between  K  and  Kp,  a  least  squares  straight  line 

was  fitted  between  K  and  Kp.  It  is  known  that  the  equation  for  the  least  squares 

4 

straight  line  fit  is  not  necessarily  reversible.  Therefore,  an  average  of  the  two 
straight  line  fits  (y  =  mx  +  c  and  x  =  m'y  +  c')  is  given  here. 

Kp  =  1.  09K  -  0.43  (1) 


4.  Spiegel,  M.R.  (1969)  Theory  and  Problems  of  Statistics,  Schaum  Publishing 
Co.  ,  New  York,  p.  229. 


In  F.q,  (l)  the  slope  is  elosc  to  1.  The  interc  ept  is  negative,  showing  that  the 
value  K  =  0  does  not  correspond  to  Kp  =  0  .  The  standard  deviation  for  the  fit  is 
0.76.  The  correlation  coefficient  of  0.  837  is  high,  showing  a  good  agreement  be¬ 
tween  K  and  Kp  for  78  to  84  percent  (0.  837^  to  0.  837)  of  the  data. 

The  magnitude  of  the  slope  and  the  high  correlation  coefficient  show  that  K  fol¬ 
lows  Kp  very  well.  Before  replacing  Kp  with  K,  one  must  look  for  the  possible 
cause  for  the  disagreement  of  about  30  percent  of  the  data  points. 

For  the  study  of  the  aurora,  the  quarter-hourly  Q  index  for  magnetic  activity 
was  suggested  by  Bartels  in  1934.  Bartels  and  Fukushima0  developed  the  scheme 
of  determining  the  Q  index.  The  Q  index  is  determined  from  stations  with  geomag¬ 
netic  latitudes  higher  than  58°.  At  present,  only  Sodankyla,  Finland,  mea  es 
the  Q  index. 

A  single -station  Q  index  is  known  to  exhibit  a  systematic  diurnal  'an  J  due 
to  the  towards-and-awav  movement  of  the  station  with  respect  to  the  aur.v,  ■  ai. 
This  is  shown  schematically  in  Figure  1,  in  which  the  Furopean  sector  is  in 

corrected  geomagnetic  coordinates.  Thu  double  ring  near  the  left  hand  be  r 

ner  shows  the  C  G  pole.  The  horizontal  line  to  the  right  is  90°  CG  longitude.  The 
geographic  North  pole  is  marked  "NP"  in  the  figure.  The  oval  rotates  around  the 
CG  pole,  while  the  station  rotates  around  the  North  pole  with  the  earth.  The  solid 
and  dashed  lines  show  the  equatorward  boundary  of  the  Feldstein  oval*1  for  Q  =  3  at 
22  and  08  UT  respectively.  The  rotation  of  the  earth  brings  Sodankyla,  Finland  ("SO" 
on  the  map)  very  close  to  the  oval  at  22  UT,  and  farthest  away  at  08  UT,  The  reli¬ 
ability  of  the  Q  index  depends  on  the  closeness  of  the  oval  to  the  station.  Thus,  the 
single-station  Q  index  from  Sodankyla  is  most  reliable  around  22  UT  and  is  least 
reliable  around  08  UT. 

For  this  reason,  mass  plots  between  Kp  (or  K)  and  Q  would  be  difficult  to  inter 
pret.  Therefore,  such  mass  plots  or  frequency  distributions  are  not  presented  in 
this  analysis. 

Differences  in  the  geographic  distribution  of  networks  used  to  determine  K  and 
Kp  might  introduce  systematic  differences  in  the  respective  magnitudes  of  K  and  Kp 
Similarly,  a  single -station  Q  index  would  generate  different  relations  between  K 
(or  Kp)  and  the  Q  index,  due  to  a  systematic  diurnal  variation  of  the  latter. 

For  operational  systems,  K  is  preferred  to  Kp  because  K  is  immediately 
available  from  the  chain  of  magnetometer  stations  maintained  by  the  Air  Force 
Global  Weather  Central.  The  magnetometer  stations  used  to  determine  the  Kp  and 


5.  Bartels,  J.  ,  and  Fukushima,  N.  (1956)  Fin  Q-Index  fur  die  Erdmagnetische 

AktivitSt  in  Viertelstundlichen  Intervallen,  Abhandl,  Akad,  Wiss.  Gottingen, 
Math-Phys.  Klasse,  Sonderheft-2. 

6.  Feldstein,  Y.I.  (1963)  On  the  morphology  of  auroral  and  magnetic  disturbances 

at  high  latitudes,  Geomagn.  Aeron.  3:183. 
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Figure  1.  Location  of  the  Q  =  3  Feldstein  Ovai  With  Hespect  to  Soda nk via  at  08 
and  22  UT 


K  indices  respectively,  are  listed  in  Table  2.  The  distribution  of  magnetometer 
networks  for  K  and  Kp  used  by  Prochaska'5  is  shown  in  Figure  2  (for  station  identi¬ 
fication  see  Table  2)  in  corrected  geomagnetic  (CG)  coordinates.  A  comparison 
of  the  distribution  of  stations  used  to  determine  K  and  Kp  respectively,  shows  that 
in  the  260°E-20°E  sector,  the  distribution  is  reasonable.  However  in  the  70°E- 
1 10°E  sector,  only  one  station  measures  K  while  s «  -en  stations  measure  Kp.  These 
differences  in  locations  and  density  of  networks  would  contribute  to  some  difference 
between  K  and  Kp.  Both  the  networks  have  no  coverage  over  a  wide  longitude  sec¬ 
tor  from  110°E  to  250°E  CG  longitude. 

To  investigate  any  systematic  differences  between  the  diurnal  variations  of  K 
and  Kp,  the  average  values  for  the  3-hour  time  intervals  and  the;r  standard 
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Table  2.  Magnetometer  Stations  for  Kp,  K  anti  Q  Indices 


Geographic 

Geomagnetic 

Symbol 

- _ 

Observatory 

Lat 

Long 

Lat 

Long 

Kp  Index 

Le 

L-erwick,  Shetland  Islands 

80  08N 

358  4 PE 

+62 .  5 

88.  6 

Lo 

l.ovo,  Sweden 

59  2  IN 

17  50E 

+  58.  1 

105.  8 

Si 

Sitka,  Alaska 

57  04  N 

224  4UE 

+60.  0 

275.  4 

Rs 

Rude  Skov,  Denmark 

55  5  IN 

12  27E 

+55.  8 

98.  5 

Es 

Eskdaiemuir,  Scotland 

55  ]9N 

358  48E 

+  58.  5 

82.  9 

Me 

Mearook,  Alberta,  Canada 

54  37  N 

246  40E 

+61.8 

301.  0 

Wn 

Wirgst.  West  Germany 

53  45N 

9  04  E 

+  54 .  5 

94.  0 

Wi 

W  tteveen,  Netherlands 

52  49N 

6  40E 

+54. 2 

9 1 .  0 

Ha 

Hartland,  Devon,  England 

51  00N 

355  3 1 E 

+54.  6 

79.  0 

Ag 

Agincourt,  Ontario,  Canada 

43  47N 

280  44 E 

+55.  0 

347.  0 

Kr 

Fredericksburg,  Virginia 

38  12  N 

282  38E 

+4P.  fi 

34  9.  9 

Am 

Amberley,  New  Zealand 

43  OPS 

172  43E 

-47.7 

2  52 .  5 

K  Index 

72489 

Houlder,  Colorado 

40  08N 

105  14 W 

+49.  0 

316.  5 

25802 

College  Observatory, 
Fairbanks,  Alaska 

84  52  N 

147  50W 

+64.  6 

2  56.  5 

71816 

Goose  Bay,  Labrador, 
Canada 

55  2  ON 

60  30 W 

+60.  5 

11.9 

727  12 

Loring  AFB,  Maine 

46  57N 

67  53 W 

+  58.  5 

1.  5 

03655 

RAF  Upper  Heyford, 

England 

51  56N 

1  15W 

+50.  7 

7  9.  1 

Q  Index 

SOD 

Sodankyla,  Finland 

67  22N 

26  38W 

+63.  4 

108.  9 

*From  Prochaska^ 


K  and  Kp  is  large  in  09-18  UT  interval.  Such  a  systematic  difference  leads  to  a 
time  dependent  relation  between  K  and  Kp. 

To  emphasize  the  large  diurnal  variation  for  a  single  station  Q  index,  the 
hourly  averages  and  average  +a  levels  for  Sodankyla  for  the  same  time  period  are 
shown  in  Figure  3  by  the  dotted  lines.  The  Q  index  scale  is  basically  the  same  as 
the  Kp  scale,  with  increments  in  steps  of  1.  The  average  curve  of  the  Q  index 
shows  a  minimum  of  1  around  07-10  UT  and  a  maximum  of  2.  3  around  22-00  UT. 
The  average  +cr  curve  also  has  a  maximum  and  a  minimum  at  the  same  times.  The 
two  curves  exhibit  a  large  diurnal  dependence  both  in  the  average  and  the  variation 
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Figure  3.  Diurnal  Variation  of  Kp,  K,  and  Q  Indices 


(a  level)  of  the  Q  index.  To  determine  the  relation  between  Q  and  Kp,  or  Q  and  K, 
the  measurements  must  be  restricted  to  the  period  22-02  UT,  when  the  statistical 
oval  is  closest  to  the  Sodankyla  station. 

For  estimating  the  extent  of  diurnal  effects,  an  average  least  squared  deviation 
straight  line  fit  (for  a  possible  failure  of  reciprocity,  see  Kef.  4)  was  obtained 
among  K,  Kp  and  Q  for  each  hour.  The  results  of  these  hourly  dependencies  are 
presented  in  Table  3.  Section  A  presents  a  correlation  between  Q  and  Kp.  Sec¬ 
tion  B  presents  the  correlation  between  Q  and  K.  Section  C  presents  the  correla¬ 
tion  between  Kp  and  K.  In  Section  C  the  correlations  are  for  every  3  hours  whereas 
in  Sections  A  and  B  these  are  hourly  values.  In  comparing  K  and  Kp  w  ith  hourly- 
values  of  Q,  it  is  assumed  that  the  3 -hour  indices  are  valid  in  each  1-hour  sub- 
interval. 

In  each  section  the  first  column  presents  the  interval  of  the  day  and  the  second 
column  presents  the  correlation  coefficient.  The  last  two  columns  present  the 
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Table  3.  Summary  of  Correlation  Analysis 


slope  and  the  intercept  of  the  empirical  least  squares  straight  line  fit.  The  last 
two  lines  list  the  highest  and  the  lowest  values  respectively  for  correlation  coef¬ 
ficient,  slope  and  the  intercept.  Table  3  (and  Figure  3)  show  a  strong  diurnal 
change  in  the  relations  among  K,  Kp  and  Q. 

The  magnitudes  of  the  correlation  coefficients  from  Table  3  are  plotted  in 
Figure  4.  In  Figure  4  the  abscissa  is  the  time  in  UT  and  th'-  ordinate  is  the  cor¬ 
relation  coefficient.  The  continuous  histograms  are  for  K  vs  Kp.  The  dashed 
histograms  are  for  Kp  vs  Q  and  the  dotted  histograms  are  for  K  vs  Q. 

For  K  vs  Kp  the  correlation  is  high  in  all  intervals,  but  dips  slightly  between 
18  and  21  UT.  The  results  show  that  K  tracks  Kp  for  most  of  the  hours,  but  this 
tracking  is  somewhat  degraded  in  the  same  18-21  UT  sector. 

The  correlation  between  Kp  and  Q  is  always  z  0.7  showing  that  these  variations 
are  also  well  correlated.  The  effect  of  the  diurnal  variation  of  Q  on  the  actual  re¬ 
lationship  will  be  considered  later. 

The  correlation  between  K  and  Q  is  z  0.68.  The  significance  of  the  variations 
is  seen  in  the  next  figure. 

Figure  5  shows  the  variation  with  UT  of  the  slopes  of  the  empirical  straight 
line  fits  from  Table  3.  The  corresponding  range  of  modulation  (highest  minus  the 
lowest  slope)  is  shown  on  the  right  hand  side.  The  modulation  is  least  for  Kp  vs  K 
and  is  highest  for  K  vs  Q. 

Before  interpreting  the  meaning  of  slope  and  its  variation,  it  should  be  remem¬ 
bered  that  K  and  Kp  have  identical  scales  ranging  from  0  to  9Q  in  steps  of  1/3.  The 
Q  scale  ranges  from  0  to  11  in  steps  of  1.  K  and  Kp  are  indices  over  an  interval  of 
3  h.  The  Q  index  has  a  much  finer  resolution  (compared  to  that  of  K  or  Kp)  of 
15  min.  A  difference  in  the  distribution  of  the  stations  contributing  to  the  respective 
indices  (Figure  2)  must  also  be  taken  into  account. 

Magnetic  disturbances  are  known  to  occur  most  often  in  the  magnetic  midnight 
sector.  Therefore  any  reference  to  time  should  be  with  respect  to  magnetic  mid¬ 
night.  However,  each  station  in  the  AFGWC  magnetic  network  has  a  different 
magnetic  local  time  for  the  same  Universal  Time.  Therefore,  it  is  convenient  to 
use  UT  in  the  following  discussions. 

The  slope  of  Kp  vs  K  is  <  1  in  the  03-09  UT  interval,  showing  that  the  response 
of  K  to  magnetic  activity  is  stronger  than  the  response  of  Kp  in  this  period.  In  me 
15-21  UT  sector,  the  slope  of  Kp  to  K  is  highest.  This  shows  that  the  response  of 
K  to  increased  magnetic  activity  has  become  weaker  in  this  period  than  the  response 
of  Kp.  The  variation  of  slope  during  a  period  of  24  h  shows  a  systematic  diurnal 
change  in  the  relationship  between  K  and  Kp.  The  03-09  UT  interval  when  K  is 
stronger  than  Kp  corresponds  to  the  magnetic  midnight  in  the  American  sector, 
where  four  of  the  five  stations  that  measure  the  K  index  are  located.  In  the  15-21 
UT  interval  the  American  sector  is  farthest  from  the  auroral  oval;  therefore  K  is 
weaker  than  Kp. 
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U.T. 

Figure  5.  Dependence  of  Multiplication  Factors  (Slope)  on  UT  for  (a)  Kp  vs  K, 
(b)  K  vs  Q,  (c)  Kp  vs  Q 


The  Q  index  also  shows  diurnal  variation,  caused  by  the  movement  of  Sodankyla 
towards  and  away  from  the  auroral  oval  as  the  earth  rotates.  Figure  1  shows  how 
the  station  is  closest  to  the  oval  near  22-00  UT  and  farthest  around  07-09  UT.  The 
Q  index  is  most  accurate  from  22  to  00  UT  and  least  useful  from  07  to  09  UT.  A 
relation  between  Q  and  Kp  is  valid  only  during  periods  when  both  Q  and  Kp  are 
accurate. 

In  Figure  5,  the  slope  of  Kp  vs  Q  reaches  its  highest  value  between  23  and 
01  UT,  the  time  period  during  which  the  Q  index  is  most  accurate.  The  slope 
reaches  its  minimum  value  between  07  and  09  UT  when  Sodankyla  is  farthest  from 
the  oval  and  the  Q  index  least  accurate.  The  overall  modulation  factor  is  2.4, 
from  this  data  base  of  39  months.  This  factor  agrees  well  with  the  previous  study 
which  showed  a  modulation  factor  of  2.  6  from  another  data  base  of  4  years. 
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Therefore,  the  slope  given  in  Table  3  for  the  22-02  UT  interval,  when  the  Q  index 
is  most  accurate,  will  be  used  in  the  empirical  relations  to  be  derived  later. 

The  variation  in  the  slope  of  K  vs  Q  is  caused  by  two  factors: 

(i)  the  diurnal  variation  in  the  response  of  Q  (measured  at  Sodankyla) 
to  auroral  oval  activity; 

(ii)  a  minor  diurnal  variation  in  the  response  of  K  with  respect  to  the 
Kp  index,  caused  by  the  derivation  of  K  from  stations  mainly  in  the 
American  sector. 

This  combination  produces  a  modulation  of  a  factor  of  3.  1  in  the  slope.  This  ex¬ 
treme  difference  in  the  response  of  K  to  Kp  and  Q  around  06-08  UT,  compared  to 
the  response  around  20-23  UT  is  due  mainly  to  the  weak  response  of  Q  to  auroral 
activity.  The  response  is  weak  because  Sodankyla  is  farthest  from  the  auroral 
oval  at  magnetic  noon. 

Figure  6  is  a  plot  of  the  extreme  values  of  K  vs  Kp  from  Tabic  3.  These  oc¬ 
cur  around  06-09  UT  and  18-21  UT.  In  the  figure,  K  is  plotted  along  the  abscissa 
and  KP  along  the  ordinate.  The  solid  line  shows  the  slope  for  the  06-09  UT  inter¬ 
val  and  the  dashed  line  for  the  18-21  UT  interval.  For  comparison,  the  dash-dot 
line  has  been  drawn  with  a  slope  of  1  (45°),  offset  from  the  origin  to  account  for 
the  fact  that  K  =  0Q  is  absent.  Comparison  of  the  lines  shows  that  K  and  Kp  are  in 
excellent  agreement  from  06  to  09  UT.  Also,  K  and  Kp  are  in  excellent  agreement 
from  18  to  21  UT  in  the  magnetic  activity  range  1.  to  2+.  For  K  >  3,  the  diiference 
between  K  and  Kp  is  successively  larger.  These  systematic  changes  are  superim¬ 
posed  on  the  local  fluctuations  in  K  that  occur  because  four  of  the  five  stations  are 
in  a  narrow  sector. 

This  information  can  be  used  to  improve  the  derivation  of  the  Q  index  from  the 
AFGWC  K  index.  The  data  from  the  period  22-02  UT,  the  period  of  highest  cor¬ 
relation  between  Kp  and  Q,  showed  (from  Table  3) 

Q  =  1.  5  Kp  -  1.7  .  (2) 

Kp  can  be  written  as  a  function  of  K  of  the  form 

Kp  =  AK  +  B  (3) 

where  the  slope  A  and  the  intercept  B  are  time  dependent  empirical  constants 
listed  in  Table  3. 

Combining  Eqs.  (2)  and  (3)  gives 

Q  =  A:  +  BXK  (4) 
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Figure  6.  Relations  Between  Kp  and  K  for  06-09  UT  and  18-21  UT  Intervals 


is  a  time  dependent  constant 

Bj  is  a  time  dependent  slope  or  multiplying  factor  of  the  index  K, 


K  is  the  magnetic  activity  index  determined  by  the  AFGWC  network. 

The  empirical  time  dependent  constants  Aj  and  computed  from  Table  3 
are  presented  in  Table  4.  The  first  column  in  Table  4  lists  the  time  interval  for 
K.  The  next  two  columns  present  the  slope  Aj  and  the  intercept  for  the  empir¬ 
ical  relation.  The  average  values  of  A^  and  Bj  are  listed  in  the  last  line  of  the 
table.  The  slope  A  j  changes  by  60  percent  between  intervals  3  and  7.  For  the 
same  intervals  the  intercept  changes  by  a  factor  of  2. 
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Table  4.  Empirical  Time  Dependent  Factors  for 
Conversion  of  K  to  Q 


Q  =  AjK  +  Bj 

where  A^and  B^are  time  dependent  terms 

Time  Interval  of  K 

Slope 

Intercept 

UT 

A1 

B1 

00-03 

1.  54 

-1.  7 

03-06 

1.43 

-1.  6 

06-09 

1.40 

-1.  8 

09-12 

1.  53 

-2.  5 

12-15 

1.  87 

-3.  3 

15-18 

2.02 

-3.  S 

18-21 

2.  03 

-3.  2 

21-00 

1.  77 

-2.2 

Average 

1.  70 

-2.  5 

Q  = 

1.  5  Kp  -  1. 7 

4.  CONCLUSIONS 


The  K  index  shows  a  skewed  distribution  with  respect  to  the  Kp  index.  At  the 
lower  end  K  is  1/3  unit  larger  than  Kp  and  at  the  upper  end  K  is  1/3  unit  smaller 
than  Kp.  In  addition  K  shows  a  diurnal  variation.  These  shortcomings  in  the  data 
are  corrected  by  deriving  empirical  relations  between  K  and  Kp. 

An  empirical  time  dependent  relation  between  Q  and  K  has  been  determined 
for  specifying  the  auroral  ionosphere  from  AFGWC  K  observations.  K  tracks  Kp 
very  well  in  the  03-09  UT  interval,  whereas  the  relationship  between  the  two  is 
somewhat  degraded  in  the  18-21  UT  interval.  Due  to  the  diurnal  motion  of  a  mag¬ 
netometer  station  with  respect  to  the  auroral  oval,  caution  is  required  in  a  direct 
use  of  the  Q-index  from  a  single  station,  like  Sodankyla.  The  overwhelming  con¬ 
sideration  is  the  time  dependence  of  K  due  to  a  limited  network  of  stations.  Statis¬ 
tically,  the  differences  between  K  and  Kp  are  negligible  for  K  =  1_  to  2+.  The  dif¬ 
ferences,  which  are  time  dependent,  become  successively  larger  for  K  >  3.  This 
time  dependence  can  be  corrected  by  use  of  empirical  factors  in  Table  4  for  a 
determination  of  the  Q  index  for  specification  of  the  auroral  ionosphere  from  K. 

The  best  longterm  approach  is  to  determine  and  correct  the  causes  (possibly 
algorithms  of  K)  to  eliminate  the  systematic  discrepancies  between  K  and  Kp. 
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Dandekar1  used  a  four  year  data  base  for  studying  the  Q-index  behavior  from 
2 

Sodankyla.  Dandekar  has  shown  that  for  ±3  hours  around  CG  midnight,  the  equa- 
torward  precipitation  boundary  determined  from  the  continuous  (diffuse)  aurora  is 
within  1°  of  the  latitude  predicted  from  the  Feldstein/Starkov^  analysis  of  Q  >  2. 
Dandekar's  results  based  on  a  larger  data  base  are  summarized  in  Appendixes  A 
and  B  for  comparison  with  results  of  the  present  study  with  a  smaller  data  base. 

Appendixes  A  and  B  consolidate  the  new  information  relating  to  magnetic  activ- 

O 

ity  indices  relevant  to  the  ERS  data  analysis. 


I'  ; 


7.  Starkov,  G.V.  (1969)  Analytical  representation  of  the  equatorial  boundary  of 

the  auroral  oval  zone,  Geomagn.  Aeron.,  (Eng.  ed),  9:614. 

8,  Gassman,  G.  J.  (1973)  Analog  Model  1972  of  the  Arctic  Ionosphere,  AFCRL- 

TR -73-0151,  AD  762  280. 


23 


References 


1.  Dandekar,  B.S.  (197  9)  Magnetic  Disturbance  Statistics  from  a  Single  Station  Q 

Index  Applied  to  an  Actual  OTH-B  Radar  Situation,  AFGl- TR -79 -02 96 . 

At)  A084  808. 

2.  Dandekar,  B.S.  (1979)  Study  of  the  Equatorward  Edge  of  the  Auroral  Oval  from 

Satellite  Observations,  AfgL-TR -79-0010,  AD  A072  997. 

3.  Prochaska,  R.D. ,  Capt. ,  Geomagnetic  Index  Calculation  and  Use  at  AFGWC, 

AFGWC/Tn-80/002.  - 

4.  Spiegel,  M.R.  (1969)  Theory  and  Problems  of  Statistics,  Schaum  Publishing 

Co.,  New  York, p.  239. 

5.  Bartels,  J.  ,  and  Fukushima,  N.  (1956)  Ein  Q-Index  fur  die  Erdmagnetische 

Aktivitat  in  V iertelstundlichen  Intervallen,  Abhandl,  Akad,  Wiss.  Gottingen, 
Math-Phys.  Klasse,  Sonderheft-2 . 

6.  Feldstein,  Y.  I.  (1963)  On  the  morphology  of  auroral  and  magnetic  disturbances 

at  high  latitudes,  Geomagn.  Aeron.  3:183. 

7.  Starkov,  G.V.  (1969)  Analytical  representation  of  the  equatorial  boundary  of 

the  auroral  oval  zone,  Geomagn.  Aeron.  ,  (Eng.  ed),  9:614. 

8.  Gassman,  G.  J.  (1973)  Analog  Model  1972  of  the  Arctic  Ionosphere,  AFCRL- 

TR -73-0151,  AD  762T50r^  - 


25 


FHSCIfilMO  1*0*  |r  K.  yj,  nim> 


Appendix  A 


Determination  of  Q  (or  Kp)  From  Photometric 
Auroral  Observations 


For  an  effective  operation  of  the  experimental  radar  system,  one  needs  iono¬ 
spheric  maps,  to  provide  an  assessment  of  the  environment  in  which  the  radio 
waves  must  propagate. 

One  of  the  important  parameters  needed  for  a  generation  of  ionospheric  maps 

of  the  surveillance  area,  is  the  magnetic  activity  index  Kp  or  Q.  The  conventional 

methods  are  unable  to  yield  values  of  Kp  or  Q  on  a  real  time  basis. 

Gassmann  suggested  the  use  of  the  equatorward  boundary  of  diffuse  aurora  as 

8 

a  real  time  indicator  for  the  magnetic  activity  index  Q.  In  the  absence  of  a  statis¬ 
tical  analysis  of  diffuse  aurora,  he  further  suggested  the  use  of  the  Starkov  equa- 
7 

tion  for  determining  the  effective  Q  index  —  Qej{  from  the  observed  equatorward 
boundary  of  the  diffuse  aurora  in  the  magnetic  midnight  meridian. 

7  0 

Starkov  fitted  an  empirical  circle  to  the  Feldstein  oval  of  discrete  auroral 
arcs.  The  resulting  equation  is  given  by: 

0  eq  =  72  -  0.  9  Q  -  5.  1  cos  (t  -  12)  (Al) 

where 

0  eq  -  is  the  equatorward  latitude  in  geomagnetic  coordinates  of  the  discrete 
aurora  (on  the  magnetic  midnight  meridian) 

Q  -  is  the  magnetic  activity  index  Q,  and 
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t  -  is  (0°  for  the  magnetic  midnight  meridian)  the  time  in  degrees 

measured  eastwards  from  the  magnetic  meridian. 

2 

Dandekar  studied  the  relationship  between  magnetic  activity  and  the  equator- 
ward  edge  of  diffuse  aurora  observed  by  wideband  spectral  photometers  aboard  the 
satellites  of  the  Defense  Meteorological  Satellite  Program.  For  the  northern  hemi¬ 
sphere,  Dandekar  derived  empirical  relations  from  a  large  data  base  of  1P72 - 1P78. 

')  eq  (diffuse  aurora)  =  7i.2  -  1.4  Kp  -2.7  cos  (t-13)  (Ac ) 


and 


■  eq  (diffuse  aurora)  «  01*.  2  -  1.0  Q  -  2.0  cos  (t-o)  (AS) 

where 

Kp  -  is  the  planetary  index  of  magnetic  activity,  and  other  terms  are  the 
same  as  in  Eq.  (Al). 

One  can  solve  Eqs.  (Al)  -  (A3)  and  determini  tin  i  iationship  between 

^discrete  arcs'  Kp'  and  Qdiffuse  aurora’  For  the  sakl  o1'  '  Onvenience  the  equator- 
ward  boundary  of  the  aurora  determined  from  the  preceding  equations  is  shown  in 
Figure  Al.  Figure  Al  consists  of  three  sections.  In  each  section  the  abscissa  is 
the  index  of  magnetic  activity  and  the  ordinate  is  the  geomagnetic  latitude  of  the 
equatorward  edge  of  the  aurora. 

The  section  at  the  top  presents  the  dependence  of  the  equatorward  edge  of  the 
diffuse  aurora  (from  DMSP  pictures)  on  the  magnetic  activity  index  Kp.  The  mid¬ 
dle  section  shows  the  dependence  of  this  boundary  on  the  magnetic  activity  index  Q. 
The  last  section  at  the  bottom  show's  the  dependence  of  the  equatorward  boundary 
of  the  discrete  auroral  arcs  (Feldstein  oval)  w  ith  Q.  The  equations  presenting  the 
respective  relationships  are  indicated  at  the  right  hand  top  corner  of  each  section. 
From  these  figures  one  can  determine  the  latitudinal  location  of  the  appropriate 
feature  in  the  magnetic  midnight  sector  if  the  level  of  the  index  of  magnetic  acti¬ 
vity  is  known. 

On  the  other  hand,  for  ionospheric  mapping  one  can  use  these  graphs  for  ob¬ 
taining  the  essential  parameter  Q  or  Kp  from  the  equatorward  edge  of  the  diffuse 
aurora  seen  in  the  DMSP  pictures. 

In  assessing  the  ionospheric  environment  for  the  radar  operation,  workers 
have  used  diffuse  auroral  boundary/discrete  auroral  boundary  and  Kp/Q  as  inputs 
in  their  computation  procedures.  One  can  determine  the  interrelations  among 
these  parameters  by  using  the  graphs  in  Figure  Al.  For  this  purpose  the  values 
of  Kp  (and  Q),  are  read  at  each  latitude  and  are  plotted  in  Figure  A2. 
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In  Figure  A2  each  curve  is  labeled  by  two  parameters.  The  beginning  parame¬ 
ter  refers  to  the  scale  of  the  ordinate  (y-axis)  and  the  end  parameter  presents  the 
scale  for  the  abscissa  (x-axis).  For  example,  the  top  curve  shows  that,  the  equa- 
torward  boundary  of  the  diffuse  aurora  at  magnetic  activity  level  Q  =  f>,  is  attained 
by  the  discrete  aurora  when  Q  reaches  (the  higher)  level  7.  This  is  consistent  with 
the  fact  that  the  discrete  aurora  usually  lies  polewards  of  the  diffuse  aurora.  As 
another  example,  the  dashed  curve  shows  that,  the  magnetic  activity  Q  =  0  is  approxi¬ 
mately  equal  to  Kp  =  i+  (l.  33).  The  difference  in  slopes  and  the  cut-off  of  the  curves 
suggest  that  one  should  be  careful  in  using  appropriate  scales  for  the  level  of  mag¬ 
netic  activity  used  in  the  generation  of  the  ionospheric  maps. 
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Appendix  B 

Relationship  Between  Kp  and  Q<jqq 

For  determining  a  relationship  be  tv..  <  K[  th«  pim  .  ta:  >  \  .!  magnetii 

activity  and  QSqd  the  magnetic  aettvit>  index  Q  fr  *t»  hi^-h  .atitu  i«  station 
Sodankyla  (Geog.  L.at.  67. 4°N,  Geog.  long.  .  •  "f  .  <  >  r  t  <  t-  •  <  ..omagnctic  I  .at. 
63.  9°N,  Corrected  Geomagnetic  i.ong.  JOB.  '«°|-  >,  tin  n  spe  ttv»  data  for  the  years 
1965,  1968,  1971  and  1974  were  used. 

A  previous  study1  of  the  Q  index  d?ta  from  S.idnnk-  la  .  Iioved  a  systematic 
diurnal  variation  due  to  the  change  in  relative  distance  between  the  oval  boundary 
and  the  station  as  a  function  of  time  of  day.  This  systematic  diurnal  variation  was 
eliminated  by  computing  empirical  detrending  factors  for  the  Sodankyla  station. 

The  Kp  index  is  available  at  intervals  of  3  hours,  whereas  the  Q  index  is  avail¬ 
able  at  intervals  of  15  minutes  .  From  these  data  a  least-squares  straight  line  fit 
between  Kp  and  QgoD  was  °btained  for  each  hour,  from  both  the  original  Q  index 
data  and  also  from  the  detrended  Q  index  data.  These  results  are  summarized  in 
Table  Bl. 

In  the  table  the  relationship  is  of  the  form  Q  =  m  Kp  +  c  where  m  and  c  are 
computed  as  a  least  squares  fit.  In  the  table  the  first  column  lists  the  year  for 
which  the  data  were  used.  The  second  column  gives  the  nature  of  data.  As  the 
diurnal  variation  shows  a  systematic  maximum  and  minimum,  the  hour  of  the  maxi¬ 
mum  or  minimum  in  the  original  data  is  listed  in  the  next  column.  The  next  three 
columns  present  the  correlation  coefficient,  the  slope  m,  and  the  offset  (cutoff)  c 
for  the  empirical  fit.  The  table  shows  that  the  original  data  exhibits  a  maximum 


Relationship  Between  Kp  a 


<♦-«  -M 

o, 

O  v 

CM 

CO 

O’  o 

o 

O-  CM 

t- 

l>  CD 

CO 

u  wn 

4J  <H  O 

CM 

CD 

m 

rf*  O 

CO 

cm  r- 

CO 

© 

CO  © 

© 

U° 

o 

O 

d 

O  f-i 

d 

d  d 

d 

o  — • 

d 

d  d 

d 

<D 

CD 

in 

C3 

o  r- 

Tf 

CO  CM 

c- 

co  m 

o 

co  r>* 

CO 

CD 

CM 

CO 

r-  02 

CD 

CO  © 

in 

O  rf 

© 

©  © 

CO 

o  £ 

CO 

H 

© 

LO  Tji 

CO 

CM 

cd  in 

CM 

m  co 

CM 

o 

r4 

o 

o  *h 

O  r4 

o 

o  ^ 

t— * 

D-  CM  S  lO  CM  ^ 

O  CM  <J  O  CM  <J 


CO  CO  3 

o  cm  <; 


5  8 

E  E 

i“"4 

< 

<< 

%  m 

to  to 

£  s 

*-H 


.  4Kp  -  0.  9  Or  Kp  =  0.  73  Q  +  0.  RG 


and  a  minimum  in  the  diurnal  variation  around  22  and  07  UT  respectively.  The 
correlation  between  Kp  and  Q  is  >  0.  (>  at  all  times.  The  slope  indicates  that  the 
relation  is  stronger  when  the  station  is  closer  to  the  oval,  and  is  weaker  when  the 
station  moves  away  from  the  oval,  u'hen  original  Q  index  data  are  used.  For  each 
year  the  ratio  of  the  slope  (m)  at  the  maximum  (22  UT)  with  respect  to  that  at  the 
minimum  (07  UT)  is  computed.  These  are  presented  in  the  last  two  columns  respec¬ 
tively.  For  each  year  the  respective  ratios  (of  in's  and  c's)  are  approximately 
equal.  The  equality  of  these  ratios  indicates  that  the  empirical  detrending  terms 
would  be  multiplying/ dividing  factors  and  could  not  be  addition,/ subtraction  type 
terms.  The  slopes  of  detrended  data  when  compared  with  the  maxima  of  original 
data  are  in  good  agreement.  The  close  agreement  shows  that  the  empirical  de¬ 
trending  factors  effectively  remove  the  systematic  diurnal  variation.  The  averages 
from  each  of  the  4  yr  are  listed  in  the  middle  section  of  the  table.  From  the  4  yr 
data  base  the  magnitude  of  the  highest  detrending  factor  is  2.  8  (for  data  around 
07  UT).  The  empirical  equations  are  presented  at  the  bottom  section  of  the  table. 
The  empirical  relations  will  enable  the  determination  of  a  Q  index  from  the  Kp 
index.  The  Q  index  data  are  needed  in  the  ionospheric  models  for  predicting  or 
estimating  the  locations  of  the  auroral  oval,  auroral  E  layer  and  the  F  layer  trough 
features. 
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